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We examine the future sensitivity of CoGeNT to the presence of dark matter streams and find 
that consideration of streams in the data may lead to differences in the interpretation of the results. 
We show the allowed particle mass and cross section for different halo parameters, assuming spin- 
independent elastic scattering. As an example, we choose a stream with the same velocity profile 
as that of the Sagittarius stream (and in the Solar neighborhood) and find that, with an exposure 
of ~ 10 kg year, the CoGeNT results can be expected to exclude the SHM-only halo in favor of an 
SHM+stream halo at the 95% (99.7%) confidence level provided the stream contributes 3% (5%) 
of the local dark matter density. The presence of a significant stream component may result in 
incorrect estimates of the particle mass and cross section unless the presence of the stream is taken 
into account. We conclude that the CoGeNT experiment is sensitive to streams and care should be 
taken to include the possibility of streams when analyzing experimental results. 



I. INTRODUCTION 

The nature of the dark matter (DM) that comprises 
the bulk of the mass in the Universe is one of the longest 
outstanding problems in all of physics. Weakly Inter- 
acting Massive Particles (WIMPs) with weak scale cross 
sections and masses in the GeV to TeV range are among 
the best motivated dark matter candidates. Recent ex- 
periments indicate hints that the DM particle may have 
been detected. However, interpretation of the data de- 
pends on our understanding of the velocities of these par- 
ticles as they pass through the detectors. The canonical 
distribution of WIMPs in the Halo of our Galaxy is the 
isothermal Maxwell-Boltzmann distribution. However, 
the real Galaxy is not likely to be quite so simple. Nu- 
merical simulations suggest that dark matter halos may 
be triaxial and anisotropic, which can lead to significant 
differences in the scattering rates and the amplitude and 
phase of the annual modulation (TJ [2] . 

Direct detection experiments are sensitive to the form 
of the local dark matter velocity distribution [3HE| , and 
the results may be significantly modified if a cold stream 
of dark matter particles exists in the solar neighborhood. 
Tidal disruption of dwarf satellites is expected to produce 
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such cold streams. The late infall of dark matter onto 
the Galaxy is also expected to result in cold streams con- 
tributing ~ a few percent to the local dark matter density 
7 - 9 . Cosmological N-body simulations from the Aquar- 
ius project pUl [IT] show the presence of tidal streams in 
the solar neighborhood at the level of ~ 1% of the local 
dark matter density, at ~ 20% probability. Larger stream 
contributions are likely if dense tidal streams pass near 
the sun's location. Streams will tend to dominate in the 
outer halo, rather than the inner halo - as has been found 
in stellar surveys thus far (for a discussion, see [T2"]). 

Recently, [T3] have studied the effect of major merg- 
ers and conclude that massive (> 10 10 Af Q ) mergers can 
lead to observable structure in the velocity distribution 
of dark matter particles. Authors [131 [J3] find that such 
a merger event is very likely. Previously |16H18j showed 
that in principle such streams should be visible in data 
from DM detectors by giving rise to a step in the en- 
ergy spectrum of the count rate: the count rate would 
be higher up to a critical energy above which the stream 
could no longer contribute to the data. Future experi- 
ments will be sensitive to a wide range of stream veloci- 
ties. In particular, as we will show, the CoGeNT exper- 
iment may have, over the next few years, the sensitivity 
to detect such streams. 

In 2010, the CoGeNT collaboration Q]5] reported an 
excess of events at low energies, which was interpreted as 
possibly due to scattering of dark matter particles with 
the target [T9H22]. More recently, CoGeNT reported [23] 
the detection of an annual modulation [231 [2H] in the 
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event rate at the 2.8a level, with a modulation ampli- 
tude of 16.6 ± 3.8%. The CRESST-II collaboration has 
also seen unexplained events that are compatible with a 
possible explanation in terms of WIMPs 26J . It is tan- 
talizing to compare these new results with the decade- 
old annual modulation seen in the DAMA data [23 EE] 
which has reached 9a confidence level. All these results 
indicate a possible ~ 10 GeV WIMP mass, yet the issue 
of compatibility between them is unclear. Some authors 
claim consistency between the different experimental re- 
sults |2"01 |2"5H52] while others do not [33H5S]. Authors 
[3"T] showed that the presence of streams could help to 
reconcile the DAMA results with other experiments. 

Yet more perplexing is the apparent discrepancy of the 
null results of XENON [351 W\ and CDMS @S1 HI] with 
both the DAMA and CoGeNT experiments; we do not 
address these null results in this paper, and restrict our 
discussion only to DAMA and CoGeNT. We further ad- 
dress the issue of compatibility between DAMA and Co- 
GeNT in this paper by extending the discussion, within 
the context of the standard Maxwellian halo, to allow 
variation in the two relevant velocities characterizing the 
distribution: the escape velocity and the dispersion. 

The focus of the paper, however, is the search for 
streams in the CoGeNT data. The low energy threshold 
(0.47 keVee) and excellent energy resolution (0.05 keVee) 
obtained by CoGeNT are key to detecting dark matter 
streams for small particle masses. The planned upgrade 
to the CoGeNT experiment (C-4) will consist of 4 detec- 
tors, of approximately 1.3 kg each [32], and is expected 
to start taking data later this year. We test the sensi- 
tivity of the CoGeNT upgrade to dark matter streams 
by performing a number of Monte Carlo simulations and 
fitting the results. We select a dark matter mass and 
scattering cross-section consistent with the current Co- 
GeNT results, assuming that the excess events currently 
seen by CoGeNT at low energies is entirely due to scatter- 
ing of dark matter particles with the target. We caution 
the reader that if a significant (and currently unknown) 
exponential background exists at low energies, our con- 
clusions may be altered. We also do not consider possible 
ways in which the CoGeNT and DAMA results could be 
made compatible with the null results of XENON [3"51l3"§] 
and CDMS [13 03]. 

Consider dark matter particles of mass m x scattering 
elastically off a nucleus of mass ton- The number of recoil 
events per unit time, per unit detector mass, and per unit 
energy, is given by the formula: 
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p x is the dark matter density at the earth's location, a p 
is the spin-independent elastic cross section for WIMP- 
proton scattering, A is the atomic mass number and 
m R,p = m x m p /(m x + m p ) is the WIMP-proton reduced 
mass. We have assumed here that the WIMP coupling 
to the proton is the same as the coupling to the neu- 
tron. F(Q) is the form factor, containing the momentum 



dependence of the cross section, and takes the form de- 
scribed in 43 -45]. f(v) is the one-dimensional speed dis- 
tribution of dark matter particles relative to the detector. 
It is this term that is sensitive to different halo models. 
^min(Q) = \J Qton/2to^, is the minimum velocity a par- 
ticle must have in order to effect a recoil at energy Q, 
where tor = m x mjq/(m x + ton) is the WIMP-nucleus 
reduced mass. 

The standard halo model (SHM henceforth) is charac- 
terized by a Maxwell-Boltzmann distribution of velocities 
in the rest frame of the Galaxy given by 
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where vq characterizes the velocity dispersion, u w h is the 
WIMP velocity relative to the halo, v csc is the escape 
velocity, and N esc is a normalization constant chosen such 
that / dvf(v) = 1. 
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FIG. 1: (a) The recoil spectrum observed by CoGeNT [23] 
and several predicted spectra compatible with the observa- 
tions. The solid (red) curve represents the best fit mass and 
cross section for i>o=220 km/s, w csc =600 km/s. The other 
panels show 3<j contours in the (m x ,<r p ) plane for CoGeNT 
(solid red) and DAMA (solid black) for (b) varying vo and 
(c) varying v csc (see text). The cross marks indicate the 
best fit points. Also shown are the exclusion limits from 
XENON10 (dashed blue), XENON100 (dashed pink), and 
from the CDMS II low energy analysis (dashed cyan). The 
XENON10 limit assumes w = 230 km/s, v CBC = 600 km/s, 
while the XENON100 limit assumes vo = 220 km/s, v CBC = 
544+®g km/s [SHIES]. The CDMS II limit assumes v = 220 
km/s, tw = 544 km/s gDHH]. 
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Fig. [TJa) shows the time averaged recoil rate for the 
low energy bins, from [23]. We have also plotted pre- 
dicted recoil spectra for WIMP masses m x = 7, 8.5, 
and 10 GeV, for the SHM assuming vq = 220 km/s and 
^csc = 600 km/s. The energy dependence of the Ger- 
manium quenching factor is obtained from the measure- 
ments reported in |46j . The solid (red) curve represents 
the best fit (m x , a p ), with a x 2 = 6.3/14 d.o.f. Panel (b) 
shows the 3a allowed contours (solid, red) for vq = 180, 
220, and 260 km/s, from left to right respectively, for an 
assumed v csc = 600 km/s. The cross marks indicate the 
parameter values that minimize the value of % 2 over the 
16 lowest energy bins. Also shown are the 3a contours 
for the DAM A results [27] [28J EI] , ignoring the possibility 
of channeling, and assuming a constant quenching factor 
for Sodium = 0.3 (the recoils off of Iodine are not signif- 
icant at these masses). We do not consider varying the 
Sodium quenching factor here, and caution the reader 
that the contours will be altered if there is a significant 
uncertainty in the Sodium quenching factor. The 90% 
exclusion limits [33(39] for XENON 100 and XENON 10 
are plotted. Panel (c) shows the variation with v esc , for 
fixed vo = 220 km/s. Shown from left to right are 3a 
contours for v csc — 600, 500, and 400 km/s. Note that 
while comparing the CoGeNT contours with the XENON 
exclusion limits, care should be taken to match the halo 
parameters. We provide CoGeNT contours for different 
values of vq and v esc , but the XENON and CDMS bounds 
are for specific values of vq and v csc . Thus the CoGeNT 
contour for v = 260 km/s is likely excluded when com- 
pared to the XENON and CDMS exclusion curves for vq 
— 260 km/s (not shown). The dark matter density at the 
sun's location was set to 0.3 GeV/cm 3 . 



II. SENSITIVITY TO STREAMS 

We now consider the possibility that in addition to 
a thermal component of dark matter, there exist dark 
matter streams, i.e. particles with small or negligible ve- 
locity dispersion. As mentioned in the Introduction, such 
streams are expected to occur due to the tidal break-up 
of small halos, the Sagittarius stream being a well known 
example. The CoGeNT experiment has also measured an 
annual modulation in the recoil rate, at the ~ 2.8er level. 
Fig. [2ja) shows the modulation (mean subtracted) with 
the expectation for theoretical models with and without 
streams. Fitting the amplitude to the SHM, we obtain 
a Xmin = 7.8/10 d.o.f, for m x = 10 GeV, a p = 0.11 fb. 
The dashed (blue) curve includes a 5% contribution from 
the Sagittarius stream for m x = 9.3 GeV, a p = 0.16 fb. 
With the Sagittarius stream included, the x min improves 
to 6.8/10 d.o.f. The pink dot-dashed curve is plotted as 
an example of an unknown stream that fits the phase 
of the CoGeNT modulation. A 5% contribution due to 
this stream improves the fit significantly, resulting in a 
Xmin — 2.7/7 d.o.f., for a small WIMP mass m x = 6 
GeV, and a very large cross section a p = 0.69 fb (for the 



unknown stream, we allow the 3 velocity components to 
vary, in addition to the mass and cross section). This 
stream has a velocity relative to the sun v sQ — (475 km/s, 
6 = 120°, 4> = 160°), and the co-ordinate system is as de- 
fined in |17j . For comparison, the Sagittarius stream's 
velocity relative to the sun [16] for vq — 220 km/s, is 
(340 km/s, 9 = 151°, (f> = 266°). We provide this ex- 
ample merely to illustrate that the presence of streams 
can have a significant effect on the phase of the annual 
modulation. More accurate results require better data 
for the annual modulation. 

Fig. [2|b) shows the 3a contours taking into account 
both the amplitude over one year, and the average recoil 
rate. The contour for the SHM is very similar to the con- 
tour obtained using the time averaged data alone (Fig. 
[T]) owing to the large error bars in the modulation data. 
The inclusion of a 5% contribution due to the Sagittarius 
stream moves the contour only slightly. The best fit val- 
ues are m x = 8.6 GeV, a p = 0.06 fb, with a % min = 14.7 
with 22 degrees of freedom (12 time bins averaged over 
energy, 12 energy bins averaged over time, and 2 fitting 
parameters). With the Sagittarius stream included (5% 
contribution), the best fit values are m x = 9.3 GeV, a p 
= 0.05 fb, with a Xmin = 14 - 9 - 

With the inclusion of a stream, the velocity distribu- 
tion f(v) is modified as: 

/(«)=£/*r(«) + (l-0/sHM(l>) (3) 

where £ is the fraction of the dark matter density con- 
tributed by the stream to the local dark matter density. 
For cold streams, f s tr{v) ~ S(v — w s tr) where v s t r is the 
stream speed. A finite velocity dispersion may be ac- 
counted for by replacing the delta function by a Gaussian. 
For definiteness, let us consider the Sagittarius stream. 
The Sagittarius dwarf galaxy is being tidally disrupted 
by the Milky Way, resulting in tidal tails. Previously, 
the possibility existed that the leading tidal tail passed 
through the local neighborhood, allowing for detection 
by direct detection experiments (for details, we refer the 
reader to [16] and references therein) . It is no longer con- 
sidered likely that the Sagittarius stream passes near the 
Earth [IS]. However, we will use the Sagittarius stream 
as a case study, to be representative of detecting a dark 
matter stream of known direction but unknown density. 
Recently, it was pointed out [JS] that the infall of the 
Sagittarius dwarf galaxy may contribute to the forma- 
tion of the spiral arms of the Milky Way. 

We study the effect of adding the Sagittarius stream 
by performing Monte Carlo (MC) simulations of future 
CoGeNT results, taking the stream to compose a fraction 
£ of the local density (which is fixed at 0.3 GeV/cm 3 ). 
We only consider the time averaged information, as the 
present data does not constrain the annual modulation 
effectively. We consider an exposure of 10 kg-year, which 
may be obtained in ~ 3 years with the CoGeNT C-4 
detector upgrade. We assume a known background that 
is modeled by a constant plus a double Gaussian, as done 
in previous works [55]. The constant and the heights of 
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FIG. 3: Monte Carlo (MC) results and fits for a £ = 0.05 
stream. The red(solid) curves show fits to an SHM+stream 
model, while the blue(dashed) curves show fits to an SHM- 
only model. The panels show the distribution of (a) mini- 
mized % 2 anc ' best-fit (b) relative stream density £, (c) dark 
matter mass m x , and (d) cross-section <r p . The vertical 
dashed line indicates the true parameters, (a) is fit by a \ 2 
distribution, while (b), (c), and (d) are fit by Gaussians. 



FIG. 2: CoGeNT modulation. Shown in (a) are the first 
12 time bins from [23], with the mean subtracted. The solid 
(red) curve is for the SHM. The dashed (blue) curve includes 
a 5% contribution from the Sagittarius stream, while the dot- 
dashed (pink) curve includes a 5% contribution from an un- 
known stream, (b) shows the 3<j contours combining mod- 
ulation information with the time averaged recoil rate, with 
and without a stream contribution, assuming no = 220 km/s 
and v e ac ~ 600 km/s. The cross marks indicate the best fit 
parameters. 



the Gaussian peaks are obtained from [53]. We include 
the known background in all our MC simulations. For 
the halo, we take vo = 220 km/s and v esc — 600 km/s 
and we take for our fiducial dark matter mass and cross- 
section m x = 10 GeV and a p — 0.05 fb. These values 
are consistent with a dark matter signal interpretation of 
the excess CoGeNT events. 

To first examine the impact of a stream on an experi- 
mental analysis, we perform 1000 MC simulations of the 
CoGeNT results for a given local stream density £ = 0.05 
(see Eq. |3|. We fit each of these simulated results to 
two types of halo models: (i) an SHM+stream model 
with variable £ and (ii) an SHM-only model (£ = 0). 
Fits are obtained by minimizing the chi-square (using the 
10 lowest energy bins in [23j ) over the mass m x , cross- 
section cr p , and, for the SHM+stream model, £. Fig. [3^ a) 
shows the minimum chi-square Xmin obtained for fits to 



the SHM+stream model (red, solid) and the SHM-only 
model (blue, dashed) . The SHM+stream model fares sig- 
nificantly better with a median Xmin °f 8.2/9 d.o.f. com- 
pared to a median Xmin °f 16.7/10 d.o.f. for the SHM- 
only model. Fig. J3^b) shows the best fit values of £ ob- 
tained for the different simulations. Figs, ^c) and (d) 
show the best fit values of m x and <r p , respectively, for 
the SHM+stream model and the SHM only model. The 
SHM+stream model gives best fit values of m x and a p 
very close to the true values. The SHM-only model on 
the other hand underestimates the mass by ps 6% and 
overestimates the cross section by ps 11%. The Sagittar- 
ius stream is clearly visible in (b) and results in erroneous 
values of m x and a p if the presence of the stream is ig- 
nored, as in the SHM-only model. 

To quantify the ability of CoGeNT to exclude the 
SHM-only halo model in favor of a halo also containing 
the Sagittarius stream, we apply a likelihood ratio test, 
which in this case is equivalent to examining the statis- 
tic Ax 2 = Xmin (SHM-only)- x mm (SHM+stream), where 
the two Xmin are the minimum chi-square obtained using 
an SHM-only halo (minimized over m x and <r p with fixed 
£ = 0) and an SHM+stream halo (minimized over m x , 
<7p, and £), respectively. The distribution of this statistic 
for an SHM-only true halo model, as determined from 
10,000 MC simulations, is shown in Fig. H (dashed blue). 
The distribution falls off rapidly for Ax^bove ~ 1, sim- 
ilar to a x 2 -distribution with 1 d.o.f. [53]. The distribu- 
tion is peaked about small Ax 2 as adding a stream com- 
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the true stream density. The exclusion level is shown for 
streams with velocity dispersions of v a — (solid red) 
and 15 km/s (dashed blue). The typical CoGeNT result 
is defined as the median A^ 2 as determined from MC 
simulations. In other words, there is a 50% chance that 
the CoGeNT results will exclude the SHM-only model 
at the given CL or better. The horizontal dashed line 
represents the 2a level. Thus, the Sagittarius stream is 
detectable at > 2a with a ~ 10 kg-year exposure with Co- 
GeNT, provided the velocity dispersion associated with 
the stream is low, and the stream contributes 3-5% of the 
local dark matter density. 



FIG. 4: The distribution of A% 2 for an SHM-only halo 
(dashed blue) and an SHM+stream halo with f = 0.03 (solid 
red) as determined from MC simulations (see the text). The 
vertical line indicates the median value of A\ 2 = 3.1 for the 
SHM+stream halo. The median value is exceeded in only 
5.3% of the SHM-only simulations. 
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FIG. 5: Confidence level at which the SHM-only model may 
be excluded when the true halo contains a stream of density 
fraction £. Exclusion levels are shown for streams with veloc- 
ity dispersions v a = (solid red) and 15 km/s (dashed blue). 



ponent to the fit is not expected to significantly improve 
the Xmin over the SHM-only fit. Also shown in Fig. [4] 
is the distribution of A\ 2 assuming the true halo also 
contains the Sagittarius stream with £ = 0.03 (solid red), 
as determined from 1000 MC simulations. Including a 
stream in the fit now allows a substantial improvement in 
the Xmin over an SHM-only fit, leading to a much broader 
A% 2 distribution. The median A% 2 is 3.1 for this halo, 
whereas only 5.3% of the simulations yield A% 2 > 3.1 
when the true halo is SHM-only. In this case, 50% of the 
time, the CoGeNT results can be expected to exclude the 
SHM-only halo in favor of an SHM+stream halo at the 
94.7% confidence level (CL). 

Fig. [5] shows the CL at which a typical CoGeNT re- 
sult can exclude the SHM-only halo as a function of £, 



III. DISCUSSION 

In this paper, we studied the ability of a future 
CoGeNT data set to detect the presence of dark matter 
streams. We performed Monte Carlo simulations of a 
halo that consists of both a thermal component, and a 
cold stream, and fitted 2 models to the simulations: (i) a 
halo model containing the stream and (ii) the SHM-only 
model. We then performed simulations of a fully thermal 
halo (i.e. SHM-only), and fitted the 2 models to the null 
simulations. We studied the Sagittarius stream as an 
example, and showed that for stream densities ^3 — 5% 
of the local dark matter density, the stream is detectable 
at the 2a level with an exposure of 10 kg year. Such 
an exposure is attainable by CoGeNT C-4 within ~ 3 
years. We set the particle mass = 10 GeV, and assumed 
knowledge of the stream velocity. Let us now briefly 
consider variations in these parameters. 

Varying the particle mass: The particle mass m x =10 
GeV provides an acceptable fit to the CoGeNT observa- 
tion, but lies at the high end of the mass range for Vq 
— 220 km/s. As the mass is lowered, we lose sensitivity 
to the stream (for vo — 220 km/s) and for m x < 8 
GeV, the Sagittarius stream becomes almost completely 
invisible as recoil events fall entirely below the energy 
threshold of 0.47 keVee. This is however, dependent on 
the assumed values of vq and v esc [50]. We have verified 
with our simulations that for m x > 8.5 GeV, we are able 
to detect the presence of the stream. Other high velocity 
streams should be visible for smaller WIMP masses. 

Varying the stream parameters: In order to test the 
importance of our knowledge of the stream parameters, 
we perform fits with a random component added to the 
stream velocity. The stream speed relative to the sun 
is chosen at random to lie between ±50 km/s from the 
true value, while the two angles that describe the stream 
arrival direction are chosen to lie between ±20 degrees of 
the true direction. This represents a small uncertainty in 
our knowledge of the stream parameters. We performed 
fits to 1000 MC simulations of the SHM+5% Sagittarius 
stream, with 25 such random velocities, and obtained 
Xmin values ranging from 8.4 to 12.1, with a median value 
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of 9.3 with 9 d.o.f. By comparison, knowledge of the true 
stream parameters yielded a Xmin 01 8.2/9 d.o.f. The 
SHM-only model resulted in a substantially worse fit (for 
£=0.05), with a x^in of 16.7/10 d.o.f. We thus conclude 
that an approximate knowledge of the stream parameters 
is still useful when analyzing data from experiments. 

We have shown that dark matter streams are poten- 
tially detectable by the future CoGeNT C-4 experiment. 
Ignoring the presence of streams may result in erroneous 
estimates of m x and a p . For sufficiently large expo- 
sures, the annual modulation provides additional infor- 
mation. The annual modulation is an excellent indicator 
of streams Hj , and should reveal the presence of the 
stream at a measured energy near the cutoff energy of the 
stream. Reconstructing the stream parameters for arbi- 
trary dark matter streams using a combination of the 
percentage modulation and the total number of recoils 
will enable us to understand the phase space distribution 
of dark matter in the solar neighborhood. In previous 
work [31 S] , we showed that streams can alter the phase 
and structure of the annual modulation. In certain en- 
ergy bins, the presence of the stream may only be ap- 
parent during part of the year, when the stream speed 
relative to the earth is largest. In these energy bins, the 
annual modulation due to a ~ few percent stream may 
result in an annual modulation of a few percent, compa- 
rable to the contribution of the entire Maxwellian halo. 



Observing the variation in the amplitude of the annual 
modulation in different energy bins provides valuable in- 
formation regarding the particle mass [HHHEl 53 . The 
phase of CoGeNT's annual modulation does not fit the 
SHM very well and in our next analysis we will investi- 
gate whether streams lead to a better fit. We plan to 
undertake this work at a later stage. 
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